LPT-ORSAY-04-63 



Brane cosmology with an anisotropic bulk 

A. Fabbri'''\ D. Langlois''''^'^ , D.A. Steer'^''^'^ and R. Zegers'^-'^''^ 

^ ! a) Dipartimento di Fisica dell'Universita di Bologna, 

O I and INFN sezione di Bologna, Via Irnerio 46, 40126 Bologna, Italy 

^ ■ and 

^ ■ Departamento de Fisica Teorica, Facultad de Ffsica, Universidad de Valencia, 

^ . Burjassot-46100, Valencia, Spain. 

^ ■ b) Institut d'Astrophysique de Paris, GReCO, FRE 2435-CNRS, 98bis boulevard Arago, 

75014 Paris, France. 

>■ ■ c) Federation de recherche APC, Universite Paris VII, 

^ . 2 place Jussieu - 75251 Paris Cedex 05, France. 

CN ■ d) Laboratoire de Physique Theorique^, Bat. 210, Universite Paris XI, 

^ . 91405 Orsay Cedex, France. 

^ : 
o . 

Oh: 

' Abstract 
> ■ 

• 1— I , In the context of brane cosmology, a scenario where our universe is a 3 + 1- 

^ I dimensional surface (the "brane") embedded in a five-dimensional spacetime (the 

' "bulk"), we study geometries for which the brane is anisotropic — more specifi- 

cally Bianchi I — though still homogeneous. We first obtain explicit vacuum bulk 
solutions with anisotropic three-dimensional spatial slices. The bulk is assumed to 
be empty but endowed with a negative cosmological constant. We then embed Z2- 
symmetric branes in the anisotropic spacetimes and discuss the constraints on the 
brane energy-momentum tensor due to the five-dimensional anisotropic geometry. 
We show that if the bulk is static, an anisotropic brane cannot support a perfect 
fiuid. However, we find that for some of our bulk solutions it is possible to embed 
a brane with a perfect fiuid though its energy density and pressure are completely 
determined by the bulk geometry. 
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1 Introduction 



Since the advent of brane models with non-compact extra dimensions, brane world cos- 
mology has been studied in depth (see [U 12 EI for reviews). The typical setup is one of 
a homogeneous and isotropic brane containing a perfect fluid, and the brane Friedmann 
equation can be obtained in (at least) two different ways. In the first approach ^ El, 
coordinates are chosen relative to the brane which is therefore at a fixed position in the 
extra dimension. The bulk, on the contrary is apparently time dependent and this time 
dependence induces, via the junction conditions, time dependence and hence cosmology 
on the brane. In the second approach jHl El, the bulk is static and the brane moves 
along the extra dimension with this motion now being responsible for the cosmology. As 
was shown in jHl, the symmetries (homogeneity and isotropy) of the brane impose that 
the static bulk is necessarily Sch-AdS. Furthermore these two approaches are completely 
equivalent via a coordinate transformation PI, and Birkhoff's theorem applies [HI. 

As in standard 4D cosmology, an obvious generalisation of such models is to consider 
homogeneous but anisotropic brane worlds [inilIIllI21linilIlllinilinilIZllIHllinil2Dlinil221- 
Recall that in 4D, it is consistent with the Einstein equations to study a universe con- 
taining a perfect fluid [23 • (We use the standard definition of a perfect fluid for which 
there are no anisotropic stresses or energy fluxes so that Tj^ = diag(— p, P, P, P) with p 
the energy density and P the isotropic pressure.) In other words, geometric anisotropy 
is consistent with matter isotropy in 4D. Thus it is possible, for example, to study in- 
flation due to a homogeneous scalar field and ask whether the universe isotropises. As 
proved by Wald [211, ^ positive cosmological constant will always isotropise the universe 
providing that the perfect fluid matter satisfies the weak and strong energy conditions. 
Observational constraints on anisotropic universes are discussed in [2S1I2H1- 

The majority^ of works on anisotropic brane world universes use the effective four- 
dimensional Einstein equations on the brane [2H1 



where '^Gab is constructed from the induced metric on the brane •jab, Sab is quadratic in 
the brane stress energy tensor Tab, and Sab is the projection on the brane of the bulk Weyl 
tensor. Although this equation looks straightforward, it must be interpreted with care 
since the Weyl term not only depends on brane quantities — either geometry or matter 
— but also on the bulk metric. Thus the above brane equations are not closed. More 
precisely, the crucial unknown quantity in this effective approach is Tr*^, the anisotropic 
component of Sab, defined in the decomposition 



Here is a given 4- velocity, hab = ■Jab + UaUb and p*, p* /3 and Tr*^ are the 'dark' energy 
density, pressure {Sab is traceless so it acts as relativistic matter), and anisotropic stress 
respectively. The dark momentum density g* vanishes for the Bianchi I branes considered 
here. The sign of p*/3 is also crucial for Wald's theorem and brane isotropisation |15j . 

^The exceptions, ^lEl, will be discussed below. 




(1.1) 




(1.2) 
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In the literature, in order to close the system of brane equations (jl-H) . various addi- 
tional conditions on Tr*^ have been proposed in a rather ad hoc way. Whether or not the 
universe isotropises, or indeed whether or not inflation can take place at all depends on 
the choice made. Similar comments hold for the isotropy of the initial singularity. 

In order to have a fully consistent picture, one cannot avoid specifying the bulk 
geometry. This has non trivial implications for the brane itself because, in brane world 
models, the Israel junction conditions — relating the extrinsic curvature to the matter 
on the brane — must be satisfied. For instance, suppose that the brane contains perfect 
fluid matter. Then the junction conditions will impose given relationships between the 
components of the extrinsic curvature: are these relations, together with the bulk Einstein 
equations, in fact compatible with an anisotropic brane? One of the initial motivations 
of this work was to try and answer that question. 

A first step was taken in |^. There the authors considered a moving brane in a static 
anisotropic bulk. Having solved the bulk Einstein equations, they then showed that the 
junction conditions induce anisotropic stresses in the matter on the brane. Hence the 
brane does not contain a perfect fluid. Furthermore, from the bulk Einstein equations 
and the junction conditions, it follows that this anisotropic stress can only vanish if the 
bulk is isotropic and hence the brane is isotropic! Thus, the conclusion is that "geometric 
anisotropy enforces, via the extrinsic curvature and the junction conditions, anisotropy of 
the matter fields" jT2j. Finally, note that the anisotropic stresses obtained on the brane 
in ^7] are fixed by the bulk metric. Thus, in general, they will not correspond to the 
intrinsic properties of physical matter. 

Another attempt was made in |16]. For the choice of bulk metric made there, it was 
shown that the brane can only be anisotropic if it contains a constant tension, p = P\ 

In this paper we go one step further in the same direction. We first start by construct- 
ing new anisotropic solutions of the 5D bulk Einstein equations. The full five- dimensional 
Einstein equations are much more difficult to solve when one assumes 3D homogeneity 
only, rather than 3D homogeneity and isotropy as in the first works on brane cosmology. 
However, with some additional ansatze, we can find explicit geometries that generalize 
the solutions previously obtained in the literature. 

We then introduce a brane in the anisotropic bulk geometries and study the constraints 
on the brane energy-momentum tensor. In particular, our aim is to see whether perfect 
fiuid matter on the brane is compatible with non-zero geometric anisotropy. We show 
that a geometrically anisotropic brane moving in a static bulk necessarily has a stress 
energy tensor which is not that of a perfect fiuid. Finally, we are able to construct a 
configuration in which matter on the brane is purely described by a perfect fiuid, though 
its energy density and pressure as a function of brane time are completely determined by 
the geometry. As we will explain, this is not so surprising. 

The paper is set up as follows. In section |21 we write down explicitly the full system 
of Einstein equations which we aim to solve. Then, in section |3[ we solve them in two 
particular cases: firstly we assume that the anisotropy depends only on time and not on 
the extra-dimension; and secondly we assume that the metric components are separable 
with respect to time and to the extra-dimension. In section |3] we embed an anisotropic 
but homogeneous brane in the anisotropic bulk solutions and discuss whether or not the 
brane matter can be a perfect fiuid. Conclusions are given in the final section. 
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2 Bulk equations 



We start from an ansatz for the bulk metric of the form 

3 



ds 



bulk 



_ 2Ao{t,w) 



,2A,(t,i«) 



(2.1) 



i=l 



where the coordinates span the three ordinary spatial dimensions and w is the coordi- 
nate of the extra dimension. Of course, one could start with a different choice of gauge, 
but we have found the assumptions g^j^ = 1 and g^fj. = convenient for integrating the 
bulk Einstein equations. 

In the bulk, endowed with a (negative) cosmological constant A = — 6/£^, the five- 
dimensional Einstein's equations simply read 



Rab = -^^gAB 



(A = 0,...,4). 



(2.2) 



Inserting the metric ansatz (j2.ip into Einstein's equations yields the following system of 
equations: 



/x=0 



i=l 
-2Ao 



1=1 1=1 

3 

l_i=0 ^=0 
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(2.3) 
(2.4) 
(2.5) 
(2.6) 



1=1 



4 = 1 



i=l 



where a dash denotes a derivative with respect to w, and a dot one with respect to t. 
Note that middle of the alphabet latin indices label the three ordinary spatial directions 
in the bulk, whilst greek indices also include the time component. In order to solve these 
equations, it is convenient to introduce the average scale factor s = (exp(^j Aj))^/^ of 
the ordinary spatial directions x*, and define 



s = e 



(2.7) 



Furthermore we denote the deviation from isotropy by a vector B with components 

Bi = Ai- A, (2.8) 

so that 

3 

J2B^ = 0. (2.9) 



1=1 



Thus B' and B quantify the spatial and temporal shear of the metric (j2.1|) . 

The Einstein equations can now be reexpressed in terms of the isotropic and anisotropic 
quantities we have just introduced. Eq. ()2.3|) becomes 



-3A - - B^ + 3AoA 
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(2.10) 



while Eq. fl2.4p can be decomposed into an "isotropic" part, obtained by averaging over 
the index i, 



-A - 3A^ + AqA 
and an "anisotropic" part that can be written as 



Be 



Finally, equations ()2.5p and ()2.6|) give respectively 



K'^ + 3A" + a;,' + 3A" + B'' = 



£2 



and 



1. 



A + (A' - a;,) A + -B • B' = 0. 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



This rewriting of the equations easily allows us to see how anisotropy modifies the usual 
equations for the bulk: setting the anisotropic parts B to zero immediately yields the 
same system of equations as in isotropic brane cosmology [HE]. 

It is thus useful to try to generalise the integration of Einstein's equations in the 
context of isotropic brane cosmology to the present anisotropic case. A combination of 
(j2.1(J|) . (j2.1H) and (j2.13|) . in fact proportional to the component Gg of the Einstein tensor, 
gives 
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£2 



(2.15) 



Another combination of the same equations, corresponding to the component of 
Einstein's tensor, gives 



A'2 ^ V ^' _ iB'2 + e-2^° i- 'A - 2A^ + ioi - ^B^^ = |. 

6 \ 6 / £2 

Following [S] let us now introduce the quantity 

Then on using ()2.6p . one can rewrite ()2.15|) . after multiplication by A', as 

+ -B'^A! + -e-2^» ( B^A' - 2A B.B' 
3 3 V 

Similarly, ()2.1(i|l yields 



Fe-^^ - -B''i - -e-2^°B2i + -A B.B' = ^A. 
3 3 3 £2 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



Note that when B = 0, one recognises the resuhs of isotropic brane cosmology in 
which case the two above equations can readily be integrated. In the anisotropic case, 
these equations are much more difficult and, in order to integrate them, we will have to 
resort to simplifying assumptions as discussed in the next section. 



3 Exact bulk solutions 

In this section we make, in turn, two particular assumptions about the metric: these will 
enable us to integrate Einstein's equation explicitly. First we assume that the shear does 
not depend on the extra-dimension, namely that B' = 0. Then we consider the situation 
in which the metric is separable, that is all the metric coefficients can be expressed as 
the product of a function of t with a function of w. 



3.1 Case B' = 

When the anisotropic parts are only time- dependent, i.e. B' = 0, then the equations 
established in the previous section simplify greatly. This situation is similar to that 
studied in ^Uj in the case of Weyl metrics. 

With B' = 0, Eq. ()2.12j) can be integrated to give 

B = A(«;)e^o-3^ (3.1) 

where A(w) is an arbitrary function of w. On substituting into (|2.18|) and (|2.19|) we 
obtain 

F'e-^^ + ^A^e-^^A' = ^A' (3.2) 

and 

Fe-'^ - h'e-'^A = li. (3.3) 
Since A is a function of w only, equation ()3.3|) can be integrated in time to yield 

F+^X'e-'^ = ^^e'^ + C{w), (3.4) 

where C{w) is an arbitrary function of w. Substitution of this ffist integral in ()3.2|) gives 
a consistency relation between the integration functions A(ty) and C{w), 

(A.A' - 2A'\^) = C (3.5) 

3 

or alternatively, using the relation A = Be~^°"'"^'^, 

e-^ {A' - A',) = 3^ (3.6) 
A {w) 

Since A and C are only ty-dependent while A and Aq in general are time-dependent, a 
simple way to satisfy this constraint is if both sides of ()3.6|) vanish. The more general 
case does not seem to yield more non-trivial solutions. 



Assuming, therefore, that C is a constant and A' = A'q, equation ()2.14|) imphes A' = 
and therefore 

A = a{t)+A{w), Ao = ao{t) + A{w). (3.7) 
This means that the metric is of the form 



+ dw'^ (3.8) 



= a\w)hf,^{t)dx^dx'' + dw\ (3.9) 

which shows that the dependence on the extra-dimension reduces to a single warping 
factor. In this sense, the anisotropy will be the same on all slices w = const, simply 
rescaled by this warping factor. 

We now solve explicitly for the warping factor a{w) = exp^(w). From (j2.5|) . or 
equivalently ()2.13|) . it satisfies the equation 

a" = ^a, (3.10) 

so that 

a'^ = ^a^ + c, (3.11) 
where c is an integration constant. Depending on the sign of c there are three solutions 

y/c£smh{w/£) c> 0, 
a{w) = { ^/^£cos^w/£) c < 0, (3.12) 
exp{w/£) c = 0, 

where we have absorbed another integration constant in a redefinition (by translation) of 
w, and w is also defined up to a sign. The other components of Einstein's equations can 
be rewritten in the form 

= (^)i?^, - (aa" + 3a'') V = V- (3-13) 

where ^^^R^u is the Ricci tensor for the four- dimensional metric h^^. This implies 

(^)i?^, = (3c) V (3-14) 

so that h^y satisfies the 4D Einstein equations with cosmological constant 6c. 

We now turn to the time dependent functions aiit) in ()3.8j) . Note that one can set 
ao = by a redefinition of the time coordinate, and this will be assumed below. For 
arbitrary c, let 

ai{t) = a{t) + (3i{t) , (3.15) 

where a{t) = ^j«j(i)/3 and '^^Piit) = in analogy with equations (j2.7ll2.9|) . Equation 
(|3.1|) can now be written in the form 

A = &^e-=^"W (3.16) 
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where the integration constants bi satisfy bi = 0. The remaining Einstein equations 
reduce to 

3a + + b^e~^" = 3c and a + 3d^ = 3c 



as given by ()3.14|1 . so that 

a2 = c + b2e-^"/6. 
On integrating it follows that, for 7^ 0, 



f: sinh(3yct) c> 0, 



and the Pi are given by ()3.16|) 



6c 

sm{3^/^t) c < 0, 
^\3t) c = 0, 



lntanh(|A/ct) c> 0, 
Intan(lv^t) c < 0, 
Int c=0. 



(3.17) 
(3.18) 



(3.19) 



(3.20) 



Here we have introduced a renormalised bi, related to bi by 

~6 



bi = 




bi, 



so that 



(3.21) 



(3.22) 



Putting together the dependence on the extra-dimension ()3.12j) and the time depen- 
dence (j3.19p and (j3.2Uj) . the full metric in the simplest case when the effective cosmological 
constant on the brane c vanishes is given by 



with 



(3.23) 



(3.24) 



We have replaced the bi by = (1 + 6i)/3, so as to recognise a warped version of the 
usual 4D Kasner metric p!Hj . 

When the effective cosmological constant is positive, c > 0, the full metric can be 
written, after some appropriate rescalings, in the form 



ds^ = smh\w/i) 



-dt^ + ^ sinh2/^(3t/£) ( tanh 



26i/3 



[dx'y 



+ dw"^. (3.25) 



Finally, the case of a negative cosmological constant is simply obtained by converting the 
hyperbolic functions into trigonometric ones. 
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3.2 Separable solution 

Let us now go back to Einstein's equations in the form (j2.3ll2.6|) and look for separable 
solutions 

A^{t,w) = a^{t)+A^iw). (3.26) 

This time we do not assume that the spatial variation of the anisotropy, B', vanishes and 
the solutions of the previous subsection, which eventually turned out to be separable, are 
a priori only a subclass of the separable solutions. 

In the equations (j2.3|) and (j2.4p . we further assume that the brackets involving time 
derivatives separately vanish.^ Leaving apart for now the mixed equation ()2.6|) the equa- 
tions involving spatial derivatives reduce to the following system: 

3 . 

■^':^ + ^',J2-^''^ = T2^ (3.27) 
E^; + E< = ^- (3-28) 

Although here we allow for a time dependence, these equations are exactly the same as 
those obtained from Einstein's equations when assuming a static ansatz, which is the 
situation studied in |17j. We can thus follow the method of 17j to integrate the spatial 
dependence of the metric. First introduce the quantity 



u{w) = exp [yAA ■ (3.29) 




Multiplying Eq. ()3.27p by u and summing over gives 

""-^^ = 0' (3-30) 

which admits the first integral, 

16 

+ =0 (3.31) 

where we have written the integration constant as — 167^/£^ for later convenience (we 
will also show that 7^ > 0). One can then integrate ()3.27|) to obtain 



where the are integration constants. 



priori, in the gauge ao = 0, one could consider the more general case where the brackets are equal 
to four different constants. However, it can then be shown that these constants are necessarily equal, 
and that if the resulting single constant does not vanish then one recovers the situation considered in 
the previous section. 
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When 7 = 0, equations ()3.31|) and ()3.32p yield 



A^{w) = w/^, (3.33) 

where the integration constants have been absorbed in a rescahng of the coordinates. 
This corresponds to the particular case c = obtained in the previous subsection. 
For arbitrary 7, the constants are constrained by Eqs. ()3.29p . ()3.28|) and ()3.27|) : 

implying \q^\ < 3/4. The second equality justifies the positivity of 7^. Equation ()3.31|) 
then yields 

u{w) = 7 sinh (4«;/£) , (3.35) 

where we have absorbed an integration constant in w. Note that replacing w by —w also 
gives a valid solution. Thus Eq. ()3.32|) gives 



1 47 

A^{w) = ^ 1^ \u{w) I + —q^v{w) + const , (3.36) 



where f' = so that 
Finally, we thus get 



v{w) = — lntanh(2|ti;|/£) . (3.37) 
47 



An 



Nlsmh^/\'i\w\/£)tanh^''^{2\w\/i) , (3.38) 
where N^^ are integration constants which must satisfy the constraint 

n K = f , (3.39) 

following from Eqs. ()3.29|) and ()3.35|) . Note that replacing w by —w in ()3.35|) would yield 
()3.38|) with —q^ in place of g^. 

As mentioned above, when the metric is assumed to be static, the full Einstein equa- 
tions reduce to the system 1)3.27113.2^ and the above results are enough to determine the 
full metric. Note that as discussed in the case {qo,qi) = (—3/4, 1/4) corresponds to 
Sch-AdSs. 

Here we assume that the metric is not static and thus the other Einstein equations 
must be integrated in order to determine the time dependence of the metric. Given 
the separable ansatz ()3.26|) . it is easy to see that the time-dependent part of the metric 
components, imposing the gauge ao = 0, must satisfy exactly the same equations as in 
section with c = 0. The solutions are given in ()3.19|) and ()3.20|) . with c = 0. However, 
there is also a further constraint coming from Eq. ()2.6|) which relates time and space 
derivatives: 

J2q^{b^ + 4:)=0. (3.40) 

i 

10 



Note that when all the g^'s are equal but non-zero, then condition ()3.40j) is incompati- 
ble with ()3.22p . which indicates that Sch-AdSs does not have a simple time dependent 
extension. In general, however, constraints (|3.22j) . (j3.34j) and (j3.4Up can all be satisfied 
simultaneously leading to time-dependent anisotropic solutions. 

To summarize, the bulk solutions we have obtained are described by a metric of the 
form 



tanh^^" (2w/£) dt^ + tanh^"^ {2w/(i) t^^' {dx')' 



+ dw^, (3.41) 



where the seven coefficients and Pi must satisfy the following five constraints 



(3.42) 



One can solve explicitly this above system of constraints and, in appendix A, we give the 
general solution for the coefficients in terms of two parameters. 

Interestingly, this metric "mixes" the five-dimensional static solution of ^Tj and the 
well-known four-dimensional Kasner solution. In this sense, we have found a much more 
sophisticated "warped" version of Kasner than in the previous subsection, because there 
are now four different warp factors, along the time and the three ordinary spatial direc- 
tions. 

As a final remark in this section, note that the above metric can be rewritten in the 
form 



ds' 



2z 



+ 



4(1-2 



2^2 



dz' 



(3.43) 



where z = ta.nh{2w / i) . In appendix B, we have used this form to compute the square 
of the Weyl tensor, which is gauge-invariant and thus useful to analyse the physical 
singularities of the metric. 



4 The brane 

So far we have obtained explicit vacuum solutions for the bulk with a negative cos- 
mo logical constant. In this section, we consider an infinitely thin brane embedded in 
anisotropic bulk geometries, and for simplicity study configurations with Z2 symmetry 
about the brane. We start by establishing the general junction conditions, and then 
discuss the possibility of embedding a brane with only a perfect fiuid as matter. 

4.1 The embedding and junction conditions 

Denoting the energy-momentum tensor of the brane matter as T^, the Israel junction 
conditions, in the case of Z2 symmetry, are given by 

= -y (n'^ - 1^5) (4.1) 
11 



where is the extrinsic curvature on one side of the brane. We use lowercase latin letters 
to denote the indices of the intrinsic coordinates on the brane. In general, the geometry 
of the brane can be defined by its embedding in the bulk space time, i.e. = X^^x"") 
where the x"" are the intrinsic brane coordinates. The extrinsic curvature is then given 
by 

= X^X^D^UB = \ [gAB{Xtd,n^ + X^daU^) + n^^^AB.c] , (4.2) 

where Da is the covariant derivative associated with the bulk metric qab^ is the unit 
vector normal to the brane, and X^ = dX^/dx"". 
If matter on the brane is a perfect fluid, 

T,'^ = diag(-p,P,P,P), (4.3) 

then the junction conditions ()4.H1 imply that the spatial components of must be equal 

Kl = Kl = Kl = -^p, (4.4) 

whilst the off diagonal components vanish. These conditions must be satisfied at the 
brane position for all times. 

Let us now introduce the bulk metric in the form 

dsl^,^ = -e'M''^)dt^ + J2 e^^^^'^^\dx'f + e^^'^'^^'^ dw\ (4.5) 

i 

where, in order to be more general, we have kept g^yj free. To obtain a anisotropic but 
homogeneous brane, we consider the embedding 

X^ = (tb(r),x,Wb(r)), (4.6) 

where the subscript "b" stands for brane. The coordinates r and x are the intrinsic brane 
coordinates. It is always possible to choose the time parameter r to be the proper time 
by imposing the condition 

e^Ao (^^y _ e^A. ^ 1. (4.7) 

As a consequence, the induced metric on the brane is simply given by 

i 

where Ai(r) = Ai{t\,{T) , W]^{t)) . 

The shear in the bulk also induces a shear in the brane geometry, which can be 
expressed as = diag(0,(T*), where 



ai ^ ^B,{tM,Mr)) = B,\^'^ + B[\^^. (4.9) 



d 

d^^ 
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The brane is isotropic when the shear vanishes. Similarly, the average scale factor on the 
brane, if, is given by 



dr 



^ + A'\^^. (4.10) 



For the embedding ()4.6|) and metric ()4.5|) the normal to the brane, n^, is 

- (e-^«+^^w;b, 0, e^°-^^tb) = - (^e-^^+^^u^b, 0, e'^'^ , (4.11) 



where a dot on and denotes a derivation with respect to the proper time (and not 
the bulk time as for the bulk quantities) and where the second equality follows from the 
condition ()4.7|1 . Here we have chosen the normal with a negative w-component, meaning 
that we keep the bulk space-time with w < Wb- For our particular bulk solutions, this 
condition will lead to p > 0. Using the definition ()4.2|) . one can then compute the spatial 
components of the extrinsic curvature tensor 

Ki, = 0, = (tj^j), Kl = ^ {n'gu + n^g'^) = n'A, + n^^A^ (4.12) 

The time/time component of the junction condition leads to the usual energy conservation 
equation on the brane, dp/dr + 3H{p + P) = 0, where H is the average Hubble parameter 
given in ()4.10|) . This is the familiar result in the case of an empty bulk and follows from 
the Gauss equation and ()4.1|1 . 

The junction conditions ()4.H1 together with the components of the extrinsic curvature 
tensor ()4.12|1 imply that the brane energy-momentum tensor is necessarily of the form 

T,- = dmg{-p,Pi,P2,P3), (4.13) 

i.e. all the off-diagonal terms are necessarily zero. However, notice that this stress tensor 
is generally anisotropic with the pressure depending a priori on the direction, unlike the 
perfect fluid form \4-^ - Let P denote the isotropic pressure and decompose 

P^=P + 7li (4.14) 

where the anisotropic stresses vTj satisfy vTj = 0. It is then instructive to decompose 
the spatial components of the junction conditions ()4.H) into an isotropic part. 



and an anisotropic part. 



X 



+ e-^Vl + e'^'^' ^'\x = -TTP, (4.15) 



^2 



+ e-^y 1 + e^^^wl B',\^ = ^vr,. (4.16) 

This last relation tells us that, in general, an anisotropic bulk geometry implies the 
existence of an anisotropic stress on the brane. 
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4.2 Anisotropic brane with a perfect fluid? 



We can now study the question of whether it is possible to find a bulk geometry and a 
brane trajectory such that (j4.16|) is satisfied with tTj = 0, 



g-Ao+A4^^ ^. ^ e-^^Jl + e^^^wl B'i\^ = 0, (4.17) 

that is, for a perfect fluid. In order to do this, we consider in turn the bulk geometries 
constructed previously, as well as the case of a static bulk geometry. 



4.2.1 A no-go condition: static bulk and a moving brane 

First consider a static anisotropic bulk, described by the metric 

rfsLik = -e^^'^^^dt^ + J2e'''''^"'\dxy + e^^^^^^dw\ (4.18) 

i 

Since we are interested in cosmology, the brane is necessarily moving in this bulk, 7^ 0. 
As the bulk is static, B = 0, and condition (|4.17p therefore reduces to 

B'l^ = 0. (4.19) 

However, from ()4.9j) . this implies that the shear on the brane also vanishes and the brane 
is isotropic. 

Thus we conclude that it is not possible to have an anisotropic moving brane con- 
taining a perfect fiuid in a static background of the form given in ()4.18p . Notice that in 
order to reach this conclusion it was not necessary to solve the bulk Einstein equations: 
a moving brane containing a perfect fiuid in a static bulk is necessarily geometrically 
isotropic, and the bulk is therefore Sch-AdSs [H]. 

Conversely, if one embeds a moving brane into such a static anisotropic background 
then the brane stress energy tensor is not that of a perfect fiuid: from the junction 
conditions, the stress energy tensor picks up a bulk-dependent anisotropic stress as in 
[I7j . In other words, the matter on the brane is fixed by the bulk geometry with the 
anisotropic stress on the brane, tTj oc B'j\-^. 



4.2.2 Bulk with B' = 

We now consider a bulk characterized by B' = 0, as in section ITTl Condition ()4.17|) then 
reduces to 



WbB 



X 



(4.20) 



so that there can only be a non-zero shear on the brane (see (jHH)) if the brane is at a 
fixed position in the extra dimension, 

Wb = constant. (4-21) 

As discussed in section 13. H when B' = the bulk Einstein equations impose that the 
bulk metric is separable (see (3.7)). This implies that condition ()4.15p reduces to the 
time independent expression 

p = A{w,)^, (4.22) 
14 



because A' depends only on w and not t. The brane can therefore contain a perfect 
fluid but it must be a constant tension given by ()4.22|) . The geometry on the brane is 
anisotropic — on choosing wi, such that a'^{wi,) = 1, 

^^Lne = -dr' + e^"^^^ E e'''^'\dxr (4.23) 

i 

where ^(t) and /9j(r) are given in ()3.19p and ()3.20j) . Finally, note that for these solutions 
the projected bulk anisotropic stress on the brane, Tr*^, vanishes since it is proportional 
to B"l , . 

4.2.3 Separable bulk solution 

We now assume that the bulk geometry is given by the metric ()3.41|) . that is ^4 = and 

= sinh^/^(4ti;/£)tanh''°(2w/£), (4.24) 
= smh^/\Aw/i)tcinh-'^°/^2w/i)t^/\ (4.25) 
: tanh^'+«o/=^ (2w;/£) t^'-i/l (4.26) 



Expression ()4.17|) summarises three equations, one for each value of i which, for this bulk 
metric, yield 



e "^o^b B', 



X 



Pi -1/3 J sinh(4M;b/^)' 



Clearly the three relations in ()4.27p are only compatible if the right-hand side is indepen- 
dent of i. Thus the coefficients Pi and must satisfy the relation 

q^ + ^ = k(p^-l) (4.28) 



where is a constant. Combining this with the constraints on pi and qi given in ()3.42j) . 
or using the results of appendix A, implies 

go = ±^ (4.29) 

and thus k = ±-\/3/2. Remarkably, therefore, the brane can support perfect fluid type 
matter for qq = As shown explicitly in appendix A, it is possible to find sets 

of coefficients {pi, g,) that satisfy all the constraints plus the additional condition go = 
±-\/3/4. All the possible solutions are expressed in terms of a single parameter. 

Of course, the brane cannot support any perfect fluid, in contrast with the isotropic 
case. The reason is that the anisotropic junction condition ()4.17|) determines the tra- 
jectory of the brane. Indeed, on substituting k = ±-\/3/2, this relation now becomes 

T^t, . ,„, ^f{h.w,). (4.30) 



^/TTwi i sinh(4u^b/£) 
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Thus, combining this with the expression for tb, one finds that the brane trajectory in 
spacetime is determined by integrating the first-order system 



dt 



dwb f 



(4.31) 



with some initial conditions t{0) = to and w{0) = Wq. Finally the energy density on the 
brane is read off from the isotropic junction condition, 



Notice that p is completely determined by the position of the brane in the bulk spacetime. 
In order to get its evolution as a function of r, one must solve for the trajectory from 
dOH and then substitute into 

As an illustration, let us consider the case go = — "\/3/4 which corresponds to a brane 
moving towards increasing values of w according to (j4.3(J|) (we implicitly assume that the 
bulk is endowed with coordinates t and w that are positive). Solving for the cosmological 
evolution on the brane, one finds that as r — oo, p converges towards the Randall- 
Sundrum tension ctrs = Q/{n'^t}. This is shown by the solid line in figure 1 which plots 
the "effective" energy density peff = p — ctrs as a function of r. This can be understood by 
the fact that the bulk geometry, at large w, resembles the warped Kasner geometry (|3.23|) . 
We have also plotted the effective pressure Peg = -P + crsj which can be computed from p 
and the average Hubble parameter via the usual energy conservation equation. Observe 
that the effective pressure is negative but converges towards zero more rapidly than the 
effective energy, so that the effective equation of state asymptotically corresponds to that 
of non-relativistic matter. 

5 Conclusion 

In this paper we have constructed explicit anisotropic brane cosmologies. We began our 
study by solving the full five-dimensional vacuum Einstein equations with negative cos- 
mological constant for a bulk metric admitting a homogeneous but anisotropic (Bianchi 
I) three-dimensional slicing. Since the general equations are too difficult, we have spe- 
cialised our analysis to two particular cases. The first case, in which the bulk anisotropy 
was assumed to be only time-dependent, lead us to bulk solutions that are warped ver- 
sions of 4D vacuum solutions of Einstein equations with a cosmological constant. These 
solutions were already known. Then we assumed separability of the metric components 
into time and extra-dimension dependent pieces. In this way, we obtained, to our knowl- 
edge, new bulk solutions that combine the 4D Kasner solution and the static 5D solutions 



We then turned to the initial motivation for this work, namely whether or not it is 
possible to embed an anisotropic brane with only a perfect fiuid as matter into a bulk 
geometry. Somewhat to our surprise, we have found that it is possible to find such a 
configuration in some of our "hybrid" bulk geometries: this is because the anisotropic 
pressure on the brane induced by the bulk time anisotropy can be compensated by a 




(4.32) 



ofM- 
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Figure 1: The effective energy density Pcs (solid line) and pressure Peff (dashed line), 
in units of a^Si ^-nd the effective equation of state Pcs/PcS (dotted line) as functions of 
the brane cosmic time r. The plot corresponds to the initial conditions tb(0) = 0.1 and 
Wb(0) = 0.2. 



corresponding term induced by the bulk spatial anisotropy when the brane is moving. 
This compensation is possible only for particular trajectories, which implies that there is 
no longer any flexibility in the choice of the perfect fluid, in contrast with the isotropic 
case. In some sense, the perfect fluid on the brane is "imposed" by the bulk geometry. 
In other words, it means that only a very particular type of perfect fluid is compatible 
with the given bulk geometry. 

This limitation is not surprising. It is simply the consequence of having less symme- 
tries. In the isotropic case, we have a high level of symmetry which implies a generalized 
Birkhoff theorem: because of the symmetries, the motion of the brane cannot perturb 
the bulk geometry, in the same way as a moving spherical shell cannot generate gravita- 
tional waves in 4D Einstein gravity. As soon as we allow for anisotropy on the brane, its 
motion in the bulk should generate a very complicated bulk. In this respect, it is a rather 
good surprise that there exists an analytical bulk solution that allows for an anisotropic 
cosmology in a purely perfect fluid brane. 
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A Appendix 

Our purpose here is to give explicitly the solutions of the constraint equations 

Y,P^ = l, $^P-=l, $^gM = 0, E^M = i' $^g.fe + l) = 0. (A.l) 

i i H IJ, i 

One can check that one can parametrize the solutions of these five constraints as follows: 



V^ = ^r,{ct>) + \, q. = \j\-\qlr,{ct> + e)-'^-l, (A.2) 



with 



and 



n{<t>) = ^ ^ cos (/) + ^ ^ sin (j), (A.3) 
o 

/ ,\ — 3 + 3 . / . ,\ 

r2[(p> = cos</)H sin0, (A.4) 

o o 

r3{(p) = — ^(cos0 + sin0), (A.5) 

All allowed sets of coefficients are thus expressed in terms of two parameters, go and the 
angle (j). 

In the last section, we considered the particular case go = ±V^/4. Substituting in 
the above results, this yields the following parametrization: 

v6 / , X 1 v2 , , „ , v3 , . ^, 

P^ = ^r,{<j)) + -, g, = _r,(0 + e±)T— , (A.7) 

with 

= 0, e_ = TT. (A.8) 



B Appendix 

In order to analyse the intrinsic properties of a metric, it is useful to construct gauge- 
invariant quantities. A particular useful quantity in the present context, where the Ricci 
tensor is already known, is the square of the Weyl tensor. 

For the metric ()3.43p . we find that the square of the bulk Weyl tensor is given by 

CabcdC^''^'' = ^ |-^^' (-2 + v^cos(30) - v^sin(30)) ^-^^"-^(l - z')t-' 



+ 



16 9 

ygo'(-27 + 56go')(l - + 8go(27 - 80go')(l - z^) + ^{23 + 2z^ + 23/) 
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-- (9 - IQqlf'^ (1 - z^) (4go(l - z^) + 9(1 + z")) (cos(3</. + 3^) - sin(3(/) + 3^)) z-\l - z^f 
9 

Ae (9 _ iQql) (1 _ _ ^og(3^) + sm(30)) ^"5/2-2go(i _ ^2)3/2^-2| _ 

Note that the spacetime becomes singular when t — > except when 30 + 7r/4 = 0, which 
corresponds to one of the pi equal to 1 while the other two vanish. 

References 

D. Langlois, Brane world cosmology: an introduction^ Prog. Theor. Phys. Suppl. 148 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 

[7 

[9 

[10 

[11 
[12 

[13 



(2003) 181, jhep-th/0209261| 

Ph. Brax and C. van de Bruck, Cosmology and brane worlds: a review, Class. Quant. 



Grav. 20 (2003) 201, |hep-th/ 0303095| 



R. Maartens, Brane world gravity, gr-qc/0312059 

P. Binetruy, C. Deffayet and D. Langlois, Non-conventional cosmology from a brane 
universe, Nucl. Phys. B 565 (2000) 269, hep-th/9905012 

P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Brane cosmological evolution 
in a bulk with cosmological constant, Phys. Lett. B 477 (2000) 285, hep- th/9910219, 

P. Kraus, Dynamics of anti-de Sitter domain walls, JHEP 9912 (1999) Oil, 
|hep-th/9910149j 

D. Ida, Brane-world cosmology, JHEP 0009 (2000) 014, 1 gr-qc /'9912062| 

P. Bowcock, C. Charmousis and R. Gregory, General brane cosmologies and their 
global spacetime structure. Class. Quant. Grav. 17 (2000) 4745, |hep -th/00 07177| 

S. Mukohyama, T. Shiromizu and K. i. Maeda, Global structure of exact cosmological 
solutions in the brane world, Phys. Rev. D 62 (2000) 024028, [Erratum-ibid. D 63 
(2001) 0299Ql], \hep-th/9912281\ 

R. Maartens, V. Sahni and T.D. Saini, Anisotropy and dissipation in brane-world 
models, Phys. Rev. D63 (2001) 063509, g r-qc/0011105, 

J.M. Aguirregabiria and R. Lazkoz, Anisotropic brane cosmologies with exponential 
potentials. Int. J. Mod. Phys. D13 (2004) 539. ,gr-qc/0304046, 

J.M. Aguirregabiria, L.P. Chimento and R. Lazkoz, Anisotropy and inflation in 
Bianchi I brane worlds. Class. Quant. Grav. 21 (2004) 823, gr-qc/0303096| 

B.C. Paul, Inflation in Bianchi models and the cosmic no hair theorem in brane 
world, Phys. Rev. D66 (2002) 124019, gr-qc/0204089 



19 



[14] T. Harko and M.K. Mak, Anisotropy in Bianchi-type cosmologies, Class. Quant. 
Grav. 21 (2004) 1489, g r-qc/0 401069 

[15] M. G. Santos, F. Vernizzi and P. G. Ferreira, Isotropy and stability of the brane, 
Phys. Rev. D64 (2001) 063506, hep-ph/0103112 . 

[16] A. V. Frolov, Kasner-AdS spacetime and anisotropic brane-world cosmology, Phys. 
Lett. B 514 (2001) 213, g r-qc/0102064 , 

[17] A. Campos, R. Maartens, D. Matravers and C. F. Sopuerta, Braneworld cosmological 
models with anisotropy, Phys. Rev. D68 (2003) 103520, he p-th/030815 8, 

J. D. Barrow and S. Hervik, Magnetic brane-worlds. Class. Quant. Grav. 19 (2002) 
155 |gr-qc/010 9084, 

[19] A. A. Coley and S. Hervik, Brane-world singularities. Class. Quant. Grav. 20 (2003) 
3061,|gr-qc/0303003, 



[20] P. Halpern, Exact solutions of five dimensional anisotropic cosmologies, Phys. Rev. 
D 66 (2002) 027503, gr-qc/0203055j 

[21] A. V. Toporensky, The shear dynamics in Bianchi I cosmological model on the brane. 
Class. Quant. Grav. 18 (2001) 2311 |gr-qc7o i03093| 

[22] M. Bruni and P. K. S. Dunsby, Singularities on the brane aren't isotropic, Phys. Rev. 
D 66 (2002) 101301, hep-th/0207189| 

[23] G.F.R. Ellis and H. van Elst, Cosmological Models, Cargese Lectures 1998, 
|gr-qc/9812046| 

[24] R. W. Wald, Asymptotic behavior of homogeneous cosmological models in the pres- 
ence of a positive cosmological constant, Phys. Rev. D28 (1982) 2118. 

[25] E. F. Bunn, P. Ferreira and J. Silk, How anisotropic is our universe?, Phys. Rev. 
Lett. 77 (1996) 2883, a stro-p h/9605123 

[26] U. S. Nilsson, C. Uggla, J. Wainwright and W. C. Lim, An almost isotropic cosmic 
microwave temperature does not imply an almost isotropic universe, Astrophys. J. 
522 (1999) LI, astro-ph/9904251 

[27] J. D. Barrow and C. G. Tsagas, The Goedel brane, Phys. Rev. D69 (2004) 064007, 
fgr-qc/0309030 

[28] T. Shiromizu, K. i. Maeda and M. Sasaki, The Einstein equations on the 3-brane 
world, Phys. Rev. D62 (2000) 024012, gr-qc/9910076 

[29] R. Maartens, Cosmological dynamics on the brane, Phys. Rev. D62 (2000) 084023, 
hep-th/0004166, 

[30] C. Charmousis and R. Gregory, Axisymmetric metrics in arbitrary dimensions. Class. 



Quant. Grav. 21 (2004) 527, gr-qc/0306069 



20 



